When a DNA Triple helix melts: An analog of the Efimov state 
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The base sequences of DNA contain the genetic code and to decode it a double helical DNA has to 
be unzipped to reveal the bases. Recent studies showed that a third strand can be used to identify 
the base sequences, not by opening the double helix but rather by forming a triple helix. It is 
predicted here that a three stranded DNA exhibits the unusual behaviour of the existence of a three 
chain bound state in the absence of any two being bound. Such a state can occur at or above the 
duplex melting point. This phenomenon is analogous to the Efimov state in three particle quantum 
mechanics. A scaling theory is used to justify the Efimov connection. Real space renormalization 
group, and exact numerical calculations are used to validate the prediction of a biological Efimov 
effect. 



I. INTRODUCTION 



Double helical DNA is common, but under certain cir- 
cumstances DNA can form a triple helix too. The 1957- 
discovery of a 3 strand DNA remained a curiouser till the 
recognition in 1987 that a third strand DNA can actu- 
ally recognize the base sequence of the double helix even 
without opening it[l|i3- Owing to its enhanced stabil- 
ity that can affect activities like gene expression, DNA 
replication and others requiring DNA opening, triple he- 
lix kindled new hopes in therapeutic applications]^. Till 
date it has been possible to make and study triple helices 
in vitro, amidst high hopes of their relevance in vivo[5|,|6[. 
At ambient temperatures, the double helix is formed 
with classical Watson-Crick base pairing while the third 
strand forms non-classical Hoogsteen or (reverse Hoog- 
steen) base pairing with one of the other two. Triple helix 
can also be formed with DNA-RNA [3 and DNA-peptide 
nucleic acid (PNA) whose uncharged peptide backbone 
helps in stabilization of the triplet structure!^, Q- On a 
completely different front, in 1970 Efimov in his studies 
on nucleons, showed an effect, now bearing his name, in 
3-body nonrelativistic quantum mechanics, viz., the pos- 
sibility of a three body bound state where none of the 
pairs are bound and the overall size of the bound state is 
much larger than the range of the pair potentials |10l - [l4| . 
One purpose of this paper is to wed these two disparate 
systems to show that an analog of the quantum Efimov 
state is the triple helix at or above the double helix melt- 
ing point. 

The origin of the Efimov effect is the scale-free quan- 
tum fluctuation near the zero-energy threshold of two 
body binding. The effective interaction that gives the 
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three-body binding does not depend on the detailed na- 
ture of the short-range interactions and is operative out- 
side its range. After many years since its discovery, it 
is now seen in systems over various length-scales ranging 
from nucleons (halo nucleus) to atomic phy sics and ultra- 
cold atoms under Feshbach resonance p^l - [l7t . A triple 
stranded DNA near its melting point is shown to be a 
unique example from the domain of classical biology and 
might provide an affordable system for studying aspects 
of the quantum Efimov physics. 

There have been many physical, chemical and bio- 
logical studies of triplex forming oligonucleotides (TFO) 
with applications in mindjlSj. We take a thermodynamic 
point of view where the long chain limit is taken to ex- 
plore the nature of phases and transitions of a triplex. 
The case of oligonucleotides can then be formulated by 
studying the finite-length effects on the transitions of the 
infinite chain system. Our focus in this paper is mainly 
on the infinite chain limit. A scaling approach is used in 
this paper to see the origin of the Efimov effect via the 
development of a ubiquitous attractive 1/r^ interaction 
over a range beyond the short range of the pair-potential. 
The short range in the context of DNA is the hydrogen 
bond length or the base pair separation. The Efimov 
interaction requires a pairwise attraction at the critical 
threshold and therefore a quantum system in d < 2 will 
not show this because in these lower dimensions a bound 
state exists for any shallow potential. However in a poly- 
meric system if a finite melting point can be induced 
then the Efimov effect would be visible. In subsequent 
sections, we consider a few simplified polymer models to 
establish the Efimov effect in DNA. A renormalization 
group (RG) approach is used to show the Efimov-like 
three chain phase in d > 2 models and a high precision 
numerical approach is used for 1 + 1 dimensional models. 
The phase diagram for the bubble-free three chain fork 
model, with a first order duplex melting, is discussed in 
the Appendix A. Some details of the RG calculation can 
be found in Appendix B. 



II. SCALING APPROACH TO THE EFIMOV 
EFFECT 

Let us consider three Gaussian polymers interacting 
with one another through a DNA base pairing type short 
range interaction. A monomer of a chain j {j ~ 1, 2, 3) 
is identified by a length variable s measured along the 
contour of the chain with rj{s) as its position vector. 
See Fig. [T^. This (and the models in the subsequent 
sections) are coarse-grained models where what we call 
a monomer in fact represents several base-pairs. The 
interaction involves monomers with same s of different 
chains, so that the Hamiltonian for strands of length TV 
(A'' — >■ oo) can be written in a dimensionless form as|19| 
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where Vki is the short-range attractive interaction repre- 
senting the base pairing between chains k and I, and Kj is 
the elastic constant or flexibility of chain j, f3 = l/fc^T, 
with T as temperature, fcs the Boltzmann constant. The 
first term on the rhs of Eq. ([T]) represent the elastic en- 
ergy or the connectivity of the polymers. The partition 
function is given by 



Z= fvTZcxp{~l3H), 



(2) 



where / DTZ stands for the summation over all config- 
urations or paths of the three chains with appropriate 
boundary conditions. We may choose rj(0) = for all 
J, i.e., the three chains are tied at the origin at one end. 
The other ends may be free. 
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FIG. 1: (a) A Gaussian chain or a random walker in three 
dimensions. An arbitrary point is specified by the contour 
length s and its position vector r(s), with the two ends at 
s — and s = N. (b) A bubble on a two chain system. 
The extent along the contour is ^y and the spatial extent is 
5x. The ladder- like regions represent the bound states with 
base-pairing between points with same s. 

By treating the partition function as a path integral, 
the imaginary time transformation s — it converts Z to a 
propagator in quantum mechanics for three particles with 
a pairwise interaction Vfe;(rfe,r/), and the masses deter- 
mined by Kj . With a short range attractive potential of 
range vq , there is a critical threshold of the potential in 
three dimensions for two particles to form a bound state, 



and for stronger potentials, there will be a finite number 
of bound states. At the critical coupling the two particle 
system has a zero-energy bound state with infinite width 
of the eigenfunction or infinite scattering length [10-[l2[. 
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FIG. 2: (a) Two noninteracting Gaussian chains 1 and 3 sep- 
arated by a distance R. Each of these can pair with a flexible 
chain 2 denoted by a thin line. The polymers are shown with 
s as an axis. In the quantum correspondence, these polymers 
are the paths with s as the time hke axis, (b) The effective 
interaction e{R) between 1 and 3. Triple chain bound state 
(Efimov effect) occurs in the region ro < -R << ^± and ex- 
tends over the whole range for ^x — > oo. 

The DNA-quantum correspondence relates the quan- 
tum critical threshold to the thermal melting of duplex 
DNA, a continuous transition in this model, with a di- 
verging length scale[20|,[2l|. This scale can be associated 
with the transient bubbles that may form for tempera- 
tures just above the melting point Tc (analogous to the 
scattering length a of the quantum problem) or the fluc- 
tuation in size and shape of the bubbles below Tc (width 
of the wave function in the quantum case). The bubbles 
can be characterized by two scales, £,± for the spatial 
extent and 
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for the length of the bubble (Fig. [T]d) where 1/z is like a 
size exponent for polymers. For the quantum problem, 
with ^_L as the scattering length a, z corresponds to the 
dynamic exponent. The diffusive or Gaussian nature of 
the free chain for our case (or the quantum problem) 
implies 



z = 2. 



(4) 



One can see this value of z from the scale invariance of 
the elastic energy term in Eq. ([T]) . This term ensures the 
connectivity of a polymer and remains invariant under a 
scale transformation, r — >■ Ar, s — >■ A^s, with Eq. (|3]). A 
similar invariance argument yields Eq. ^ for the nonrel- 
ativistic free particle Schrodinger equation. 

Suppose we have two noninteracting chains I and 3, 
both interacting with another one, chain 2. For siinplic- 
ity, though not essential, we may take I and 3 as rela- 
tively stiffer compared to 2, and they are at a distance 
R apart (see Fig. ^ . Choose the pairwise interaction be- 
tween 1-2 and 2-3 or the temperature close to the critical 
threshold. In that case, if £^± > R, then inbetwecn two 



contacts with chain 1, chain 2 is expected to meet chain 
3. The interaction is attractive. This exchange over a 
large length scale induces an attraction between 1 and 3 
if chain 2 is integrated out from the partition function of 
Eq. ([2]). The effective interaction e{R) is the change in 
free energy AF/N because of the presence of a scale R 
so that it can be written as 
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where the first factor is the number of blobs and J^{x) 
is a scaling function. For ^j^ ^ oo, Eq. (O should be 
independent of it requiring J^{x) ^ x~' , as a; — ^ 0, so 
that in this hmit, by Eq. ^, 

e(-R) = — ^ ^-^ = ";^' (^ = a constant). (6) 

We see the emergence of a "universal" 1/R^ interaction 
for a region rg < i? <C ^± , and this attractive long range 
interaction can produce a bound-state of size bigger than 
rg. To be noted here, that above the duplex melting 
point, there will be large bubbles, thanks to critical fluc- 
tuations, and so the attractive long range interaction per- 
sists even above the melting point. This gives the possi- 
bility of a three chain bound state where none are bound 
in the two particle potential well. 

In the quantum language the large fluctuation is the 
resonance, and the Efimov effect is due to the attraction 
produced by the multiple scattering of a light particle 
off the two heavier ones. This is one aspect of the Efi- 
mov effect. In an exact study using a separable potential 
(under the Born-Oppenhcimer approximation) the effec- 
tive interaction between 1 and 3 by the hopping of 2 has 
been calculated [1^. The result of Ref. [ij] can be recast 
in our scaling language (note that ^_l is the scattering 
length a) as 



e{R) = -^ irf(i?), 
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corroborating the scaling hypothesis of Eq. ([5]) . One sees 
a cross-over from —XjR? for i? <C 1 to the Yuakawa form 
—e~^/R for R ~ 0(1). The scale free interaction per- 
meates the whole region for ^^ ^- oo at the critical or 
threshold value of the pair interaction. Fig.[2]3 shows the 
nature of the effective interaction in presence of a scale 
^j^ . The Efimov interaction is beyond the hydrogen bond 
length of duplex DNA. 

The Efimov interaction allows a three chain bound 
state close to, or just above the duplex melting point. 
The size of the bound state will necessarily be much 
larger than the hydrogen bond strength. The fact that 
one sees Efimov states, though finite in number, for large 
but finite ^x (Eq. (|S])), indicates that the triple helix 



DNA may also show an Efimov analog bound state even 
above the duplex melting temperature. This also means 
that the melting temperature is higher for a triplex than 
for a duplex. The minimal model in Eq. ([l} has a con- 
tinuous transition. If other interactions not in Eq. ([1]) 
drive the transition first order, the length scale ^j_ will 
be non-diverging. However, if the transition is weakly 
first order, ^x raay be large enough to accommodate the 
intermediate scale- free region shown in Fig.[2lD allowing a 
bound state. Hence in a weak first order transition with 
large ^j_ one would still see the Efimov-like bound state. 
A case of a strong first-order transition is dealt with in 
Appendix A where it is shown that the effect is absent if 
the bubbles are fully suppressed. 

The second aspect of the Efimov effect is that at the 
critical point, there is an infinite sequence of bound 
states in a geometric form En = Eo{e~^'^^'^°)" where sg 
is system specific number (e.g. for 3 identical bosons. 
So — 1.00624). The energy scale Eq is also system- 
parameter dependent. The above analysis is done with 
the ground-state dominance assumption that is justifi- 
able for N ^ oo. The polymer partition function of 
Eq. ([2]) for finite N is given by the Efimov energies as 

Zr^Ci expi^EnN) + C2 exp{-EiN) + ..., (9) 

where the temperature has been absorbed in the "ener- 
gies" . The terms beyond the ground state can be ignored 
if 
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This is the length requirement for the triple DNA to show 
the Efimov-like state. 

The prediction of a large girth triplex close to (and 
even above) the duplex DNA melting point awaits ex- 
perimental tests. 



III. EFIMOV-LIKE PHASE IN rf > 2 

In order to justify the prediction of the scaling the- 
ory in a systematic way, we adopt a renormalization 
group [22[ approach for d > 2. In the RG approach, 
the effects of interaction is probed by summing over the 
configurations at a smaller scale (in the partition func- 
tion) and redefining the effective interaction on a larger 
scale. For a bound state, we should see an effective in- 
teraction among the chains, irrespective of the scale of 
coarse-graining. In contrast, for an unbound state, lo- 
cally bound monomers lose their importance as we sum 
over configurations and therefore the effective interaction 
would vanish as the probing lengthscale increases. These 
effects are expressed by the RG flow equations or recur- 
sion relations, as flows of the interactions with length 
scale. A two body bound state should therefore be pos- 
sible if the two body interaction does not vanish. In the 
same spirit, a three body bound state would occur if a 



three-body interaction becomes operative, even if there is 
none to start with. We express these RG relations in an 
exact form on specially constructed hierarchical lattices 
with discrete scaling symmetry and tunable dimension- 
ality. 



or absence of three body interaction and the unstable 
fixed point b^ — 1 describes the critical state produced by 
a pure three body interaction. The flow going to infinity 
is indicative of the three-chain bound state, formed by 
the three body force, a case not of interest here. 




n=0 



n=l 



n=2 



FIG. 3: Recursive construction of tlie hierarchical lattice used 
for the real space renornialization group. At every stage, each 
line is replaced by a diamond of 2fe lines. 

Let us consider three directed polymers, which are do- 
ing random walk from bottom to top on a hierarchi- 
cal lattice, constructed recursively with a motif of 2b 
bonds, as shown in Fig. [31 The branching factor b deter- 
mines the effective dimension((i) of the hierarchical lat- 
tice as d = Y^- There are attractive potentials — e^ and 
—^ijki^ijjSijk > 0) if a single bond is shared by two and 
three polymers respectively. Although £123 = 0, still it 
will be needed for the RG transformation to probe the 
three-body bound state and it is generated by renornial- 
ization. 

The configurations of the two chain system can be clas- 
sified as two independent chains or inherently two chain 
configurations. The corresponding recursion relation for 
the two chain Boltzmann factors yij = exp(/3eij) is given 
bym 



5-1 



y^ 



yt 



(11) 



Siinilarly the three chain configurations can be classi- 
fied as (i) three independent chains, (ii) a combination of 
a double and a single chain, or (iii) inherently three chain 
configurations, i.e. three chains sharing the same bond. 
By a decimation of the 26-bond motif, the recursion rela- 
tion for the three chain Boltzmann factor w = exp(/3ei23) 
is given by 



{b~i){b-2) + {b-i)Y,yf, + ^'l[y^i' 
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i<3 



(12) 



where y',w' denote the renormalized values. 

In the following discussion, we choose the initial value 
w = I. The three fixed points of yij correspond to (1) 
y* — 00 (zero temperature, pure bound state), (2) y* = 1 
(infinite temperature, no two body interaction), and (3) 
y* = (6—1) (two chain unstable critical point). 

In case there is no pairwise bound state or no pair in- 
teraction, w has two fixed points 1 and &^ — 1. The stable 
fixed point 1 describes the high temperature fixed point 



(a) 




FIG. 4: (a)RG Flow-diagram in the y-w plane for the symmet- 
ric case j/12 = 2/23 = yis, all starting with w = 1. Here 6 = 4. 
The flow of w goes to cxj if the starting y > yo = 2.32402..., 
otherwise to 1 (high temperature flxed point). The trajecto- 
ries with starting y < yo end sd, w = y — 1. (b) Average 
energy per monomer vs. temperature from direct computa- 
tion (chain length=2^'''). For two chains (marked 2) average 
energy undergoes a continuous transition ai y = yc while the 
average energy for three chains (marked 3) shows a jump at 
y — yo- The region from yo to j/c is the Efimov-like three 
chain bound state. 



In case all pairs are in the critical state so far as the two 
body interaction is concerned (y*- = 6—1), there is no real 
fixed point for w. The renormalization flow takes w to 
infinity, as shown in Fig. |4^. The three chains then form 
a bound state at the two-body critical point. For temper- 
atures above the duplex melting, i.e., with initial values 
y = yi2 = 2/23 = J/31 < 6 — 1, the triplex will be in the 
denatured state if the flow goes io y — l^w — 1. E.g., for 
6 = 4, the three chain melting is at y = 2.32402... which is 
less than the duplex melting point ^c = 3. A further con- 
flrmation of this triplex melting comes from an exact nu- 
merical calculation of the average energy by iterating the 
partition functions and their derivatives for large lattices. 
This is shown in Fig. St). With y = yi2 = j/23 = 2/31, as in 
Fig. 2^, the two strand system melts through a second or- 
der transition &t y = y^ (energy going continuously to 0) 
whereas the three strand system undergoes a first order 
transition at a temperature y = yo < yc (energy showing 
a discontinuity - see Appendix B for details). The region 
between y ^ yo to y ~ yc, is the region of a triple strand 
bound state when there should not be any duplex. 

The phase diagram in the plane of y^^ vs yi2 with 
2/12 = 2/23, w = 1 is shown in Fig. [5] For y^^ = 0, chains 
1 and 3 are bound for ever and therefore chain 2 melts 
off at yi2 = y/b — 1. This point is indicated by a star in 
Fig. [5] Within the triangular shaded region bounded by 
y^^ — 1/(6— 1), y^2 — l/(^ ^ 1), a-iid the curved line 
separating the unbound state, we have a triplex phase 
without pairing of any two - the desired Eflmov effect. 
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FIG. 5: Phase Diagram: y^^ 
The duplex melting point at y, 



vs. 



yi2^ (2/12 =2/23), for 6 = 4. 
b — 1 is indicated by 
the horizontal and vertical lines. Three chains are bound in 
the shaded region with the thick curve representing the three 
chain bound-unbound transition. Above the horizontal line 
at 2/31 = & — 1 in the shaded region, a triplex state exists even 
though no two should be bound. The bound states in other 
regions are as indicated. The star at y^2 ~ l/V^ ^ 1 is the 
melting of chain 2 and composite 1,3. 



IV. NUMERICAL SIMULATION IN d = 1 + 1 

After arguing for the existence of an Efimov bound 
state for triple DNA helices through the DNA- quantum 
correspondence and the scaling argument, and after ob- 
serving the effect on a hierarchical diamond lattice, in 
this section we produce a clear numerical evidence that 
the effect is present in an euclidean lattice even in 1 + 1 
dimensions. 

Directed polymers in such a dimensionality are 
amenable to exact solutions in the thermodynamic limit 
of infinite chain length or for extremely accurate numer- 
ical simulations since their interchain contacts are guar- 
anteed to occur between monomers (see fig. |6l) with the 
same index, as for DNA. For such reasons, they played an 
important role in clarifying melting, cold unzipping and 
overstretching properties of duplex DNA |23l - l25l | , and for 
the same reasons they turn out to be the most conve- 
nient models in which to test numerically the existence 
of the Efimov effect that on higher dimensions might be 
obscured by the noise of numerical simulations. 

Here, we consider a duplex DNA formed by directed 
polymers which can cross each other (see Fig. ^ . Such 
models do not show any melting transition at finite tem- 
peratures in 1 -|- 1 dimensions. In the quantum mechanics 
analogy, the dimension along which the polymers are di- 
rected plays the role of time: then it is well known that 
any short-range potential, no matter how weak, will pro- 
duce a bound state for d < 2. 

Nevertheless, keeping the notation y = exp(/3e) for the 
Boltzmann factor associated with base pair interaction, if 
we introduce a fugacity a for each bubble formed in the 
model between the two DNA strands, it is possible to 
demonstrate that the melting transition 2/c(c) decreases 
from yc{l) = 1 (i.e. (T = 00)) for ct = 1 to 2/^(0) = 2 
for a = (i.e. for the fork model in which bubbles are 
suppressed. See Appendix A for the phase diagram of 





FIG. 6: a) Graphical representation of directed two strand 
DNA in l-f-l dimensions. The two strands start from the 
same origin and are directed along the direction (1,1). The 
walks can cross each other. Each bubble opening is weighted 
with a factor a, each interaction (overlap of the two chains) 
with a weight y. The weight of this conformation is there- 
fore cT^y''. Notice that each interaction corresponds to sites 
with the same index, b) The three chain models. A factor a is 
given to each bubble opening among all possible pairs of poly- 
mers. For model A the weight of this configuration is a^^y^^ 
(see text) for model B, where the interaction is only between 
the red and the blue chain, and the red and the brown chain, 
the weight is a^^y^. The black arrows, in both figures, indi- 
cate the end to end distance rjv used to estimate the melting 
transition. 



In order to generalize the model for three chains (see 
fig. [U we associate the fugacity a to each bubble opening 
between all possible pairs of chains. As regards the base 
pair interactions, we consider two options: 

• In Model A we assign a Boltzmann factor y for 
all two chain interactions, but in the case when all 
three chains come together we assign an interaction 
factor y^ instead of y^. 

• In Model B, we assign a Boltzmann factor y for 
each interaction between chain 1 and chain 2, and 
between chain 2 and chain 3, but we do not consider 
any interaction between chain 1 and chain 3. 

In both cases, the melting temperatures of the three and 
two chain systems coincide at cr = 0, when bubbles are 
forbidden. The existence for a given cr of a three chain 
bound state at a y < yc(o'), for < cr < 1, should be a 
definitive proof of the Efimov effect for models of triplex 
DNA. 

It is known [20j that in 1+1 dimensions, if a poly- 
mer is confined between two lines or two polymers, there 
is a repulsive entropic force of a similar 1/?'^ potential. 
In our model B, chain 2 is the only one that can medi- 
ate interaction between 1 and 3, because the latter two 
chains do not interact, and for chain 2 to do so, it gets 
effectively, though not strictly, confined between chains 
1 and 3. In such a situation the steric repulsion and the 
induced attractive interaction between 1 and 3 tend to 
cancel each other [2y|. This weakens the possibility of the 
Efimov effect in model B but not in the case of model A. 



The nonexistence of the Efimov effect for model B would 
therefore be a further proof that our analysis is correctly 
taking care of all possible interactions of the models. 

The three chain models cannot be solved exactly but 
they can be numerically studied with impressive preci- 
sion through transfer matrix techniques. We worked at 
CT = 0.5. For this value the melting transitions occurs 
at y„i = j/c(0.5) = 4/3. The existence of the melting 
transition can be seen by looking at the rescaled average 
distance S^n = ff^ /{N'^/'^) between the free extrema of 
an arbitrary pair of the three chains at a length N . By 
using standard finite size scaling techniques, in order to 
pinpoint the transition temperature, we computed the in- 
tersection between these curves. The resulting numbers 
are presented in Fig. [T] In the figure, besides the data 
for model A and model B, we also show the data for the 
exactly solvable two chains model in order to evaluate 
the degree of convergence of the simulations. 
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FIG. 7: For models A and B described in the text and for 
a two strand model, we estimate the denaturation transition 
by looking, in the three cases, at the crossing yt of the curves 
^N{y) and ^jv+ioo, where A'' is ranging from 1500 to 3900 in 
steps of 100. These crossing points are plotted as a function of 
1/(A'' -I- 50). The case a = 0.5 is considered. The horizontal 
continuous line represents the exact melting value yt = 4/3 of 
the duplex model. Linear extrapolation curves are also shown. 
Their intercepts with the y-axis are 1.3333 ±0.0001 for model 
B, 1.3333 ± 0.0001 for the 2 chains, and 1.3326 ± 0.0001 for 
model A. These results definitely show that model B has the 
same melting transition as the duplex model, whereas model 
A melts at a higher temperature (Efimov effect) 

Beyond any possible numerical uncertainty, model B as 
we argued already does not show any Efimov effect while 
model A exhibits the Efimov effect. Although narrow, 
there exists a temperature interval for model A in which 
three chains are bound whereas two are not. 



V. CONCLUSION 

We have presented a scaling argument to show the pos- 
sibility of a three strand DNA bound state at conditions 
where a duplex DNA would be in the denatured state. 
This is a biological analog of the nuclear or cold atom 



Efimov effect. The scaling argument is further confirmed 
by renormalization group and exact numerical results on 
various model systems in different dimensions. 

We end with a few comments. To mimic reality, the 
minimal model considered here may be supplemented 
with additional terms like excluded volume effects whose 
influence on the phenomenon need separate analysis. 
However if experiments can be done in thcta conditions 
for the strands, then the excluded volume effects can be 
ignored or minimized, and our results would be applica- 
ble. It may be noted that enzymatic activities are hy- 
pothesized to involve local covalent or hydrogen bonds 
and contacts. For DNA, one requires denaturation of 
strands (melting in some form) locally. In these lock and 
key mechanisms, the issue of nonspecific long range bonds 
is generally not considered. Therefore, the existence of a 
bound state involving two otherwise denatured strands 
of DNA due to the presence of a third strand, with 
overall separation much larger than the hydrogen bond 
length would have important implications in biological 
processes. Nonetheless, we anticipate new experiments 
to look for signatures of our proposed Efimov-DNA. 



Appendix A: Phase diagram in absence of loops 
(fork model) 

In order to get the thermodynamic phase diagram, we 
consider a simplified model of three chains without bub- 
bles by generalizing the double chain Y-model studied 
in the context of melting and unzipping transition of a 
duplexlHlll. 

Consider 3 chains on a standard square lattice, 
stretched along one diagonal direction (directed poly- 
mers). The polymers interact with a pairwise contact 
interactions — e^j between chains i and j. The monomers 
are indexed from one end s = 0. By construction, the in- 
terchain interaction is between monomers with the same 
index, as is required for DNA. There is an additional 
constraint that an unbound configuration of two chains 
cannot be followed by a bound stretch along the increas- 
ing s direction. This avoids loops on the chains but allows 
opening of the DNA. This is the Y-model. For this model, 
every open chain has an entropy fc^ In fj, per monomer 
(/i = 2 steps per bond) and the same entropy for any 
bound polymer (duplex or triplex). With this entropy 
and the additive energy, the free energies per monomer 
for various possible phases can be written down. These 
arc (with kg = 1) 

/b = -(ei2 + e23 + e3i)-rin/i, (allbound)(Ala) 

/u = -3Tln/i, (all unbound), (Alb) 

fij,k = ^^ij — 2nn/i, {ij bound, k free), (Ale) 

The two chain melting transition is at Tc = ejfe/(ln/i). 

The phase transition lines in this Y-model are all first 

order lines whose slopes can be determined by a Clausius- 

Clapeyron argument. Let us take a line of coexistence 



(a) Tln2>c>0 



(b) c=Tln2 




bound state is purely from the pairwise binding and in 
that sense it is not a true representative of the effect j27| . 
The figures show that at or above the mehing tempera- 
ture of the 1-3 pair or if chains 1 and 3 are noninteracting 
(c = 0, Fig. [S};), the three chain bound phase may melt 
via a duplex or directly to the unbound state. Direct 
melting without an intervening duplex phase is not pos- 
sible if c< 0, (Fig.IHJi). 

Experimental phase diagrams for triple helix are gener- 
ally done with temperature and salt concentration as the 
two variables. For a given concentration, the variation of 
the temperature would follow a particular trajectory in 
the three dimensional thermodynamic phase space and 
the possible phases one would see are determined by the 
intersection of that curve with the phase boundaries. 



Tln2+lcl 



FIG. 8: Phase diagrams for 3 chain Y-model. Here c = esi. 
c < corresponds to repulsive interaction. "U" represents 
the unbound state. 



between two types of phases A and B in a phase diagram 
of "intensive" variables Xi and X2. Let the conjugate 
extensive variables be pi = dF/dXi, F being the appro- 
priate free energy. At a given point on the coexistence 
curve Pi would show a discontinuity, taking values pf and 
pf in the two phases. Then the slope of the line is given 
by 



0X2 
dXx 



Pt 



P? 



Pi 



pT 



(A2) 



In our case, X\ = /3ei2, X2 = /3e23 and so the derivative 
of the free energy gives the associated number of con- 
tacts. Consequently, the slope is related to Ani2/An23. 
Now in this all-or-nonc Y-model, these differences in the 
fraction of contacts in the two phases are either 1 or zero. 
Therefore, the slopes can in general be or 00. The spe- 
cial case is line Xi = X2 ■ By symmetry, here on this line 
Ani2 = ±Ari23 and the slope should be ±1. 

If one chain melts away, i.e., a triplex breaks into a du- 
plex and a free chain, then Tb^i2,Tb^23 < T^^u- These 
conditions follow from the stability of the phases as given 
by the free-energies. We therefore have two inequalities: 
for the {12,3}-state to occur, £12 > £23 + £31 and for 
state {23,1} to occur, £23 > £12 + £31- For easy dis- 
play, we represent the three dimensional phase diagram 
in slices of the £i2-£23 plane for given values of c = £13 
and T. Fig. [5] shows four possible situations for different 
values of £31 and T. In these figures, c = Tln2 is special 
because of melting of the 1-3 pair. There is a region in 
Fig. [8^, where one sees a three chain bound state in a 
range of temperatures where none of the pairs would be 
bound. Despite the similarity with the Efimov effect, the 



Appendix B: Evidence of a first order transition for 
triplex 

We present a numerical evidence that there is a dis- 
continuity in the three chain average energy as shown in 
Fig.H). 

The exact value of En at the nth generation is com- 
puted by using MATHEMATICA with infinite precision 
by iterating the two and the three chain partition func- 
tions and their derivatives. The length of the polymers 
at the nth generation is i„ ~ 2" so that the thermo- 
dynamic limit of energy per monomer E„/Ln can be ob- 
tained by extrapolation to 1/n — > 0. If C„, Z„ and Qn are 
the nth generation partition functions for single, double 
and triple chain systems, then these obey the recursion 
relations [22I 



Cn 


- bCLi, 




(Bl) 


Zn 


= &(&-l)Cf,_i+6Z^i, 




(B2) 


Qn 


= b{b-l){b-2)Ct^, 








+6(6-l)C,ti5]Zfe-)L 


1 + &QL 


-i,(B3) 




i<3 







where the arguments of Z„_i in Eq. (|B3p refer to the two 
chains involved. The initial conditions are 



Co = l,Zo = y,( 



y 



The relations for the derivatives can be derived from 
Eqs. (lBT])- (lB3|) . 

Fig. [5] shows the extrapolation in the range oi y = 
2.323 to 2.327 which brackets the transition in the 
range (2.324,2.325). The discontinuity survives even 
on a finer scale in Fig. ^jp, which give yo in the range 
(2.32402, 2.32403) consistent with the RG result of Fig. g] 
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FIG. 9: Plots of En/ Ln vs 1/n, En being the three chain average energy, n the generation number. We went upto n = 25 for 
which the length of each polymer is Ln = 2'^^. In (a) we show for y — 2.322 + 0.001 *n,n =0-5, while a finer grid result is 
shown in (b) with y = 2.32400 + 0.00001 * n,n =1-4. The lines show the extrapolations to n -> oo. The discontinuity at the 
transition is visible. 
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